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O Abstract 

^ We discuss the size of isospin violating effects of both electromagnetic (e 7^ 0) and 

• T-j strong (m„ 7^ m^) origin in ^^2(7) decays. Our results are relevant for the extraction of 

r> from the measured ratio of the K^j^2{^) and 7r^2(7) decay widths combined with 

Fk/Ft^ obtained in lattice calculations in the isospin limit. We point out that strong 
isospin-breaking corrections, neglected in previous studies, enter at the same level as 
the electromagnetic effects. The updated value for the ratio of CKM elements reads 
\Vus/Vud\ = 0.2316(12). 
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1 Introduction 

The purely leptonic decays of light pseudoscalar mesons P — > iu^ (P^ = n^,K^) are 
theoretically extremely clean. Neglecting electroweak corrections, the hadronic physics is 
encoded in one single parameter, the meson decay constant Fp. Using chiral perturbation 
theory and lattice QCD, P — )■ iui decay rates can be predicted with high accuracy and 
provide non-trivial tests of the standard model imposing stringent constraints on new 
physics. 

With the advent of high precision calculations of the ratio of decay constants Fx / F^^ in 
lattice QCD, it was realized [1] that the combination of Fk/F^^ and the ratio of measured 
decay widths r_ft'^2{7) /-'^ 1-^2(7) ^^^^^ ^ clean and competitive determination of the ratio 
of Cabibbo-Kobayashi-Maskawa (CKM) [21 [3] matrix elements. Implementing 
the programme proposed in Ref. [1] requires theoretical understanding of: (i) the long- 
distance electromagnetic (EM) corrections to P — )■ iuf, (ii) the strong isospin breaking 
correction that connects Fk/F^^ (calculated within lattice QCD in the isospin limit m„ = 
rrid) to Fk± I Pn-i that enters the ratio of K^2{-y) to vr^2(7) decay widths. 

The EM corrections have been known for some time [U El El El • The analysis of Ref. [T] 
relied on the model-dependent treatment of Ref. [7j for the long-distance contribution 
to these corrections. The long-distance EM corrections to P — > ^2/^(7) have also been 
analyzed in the model-independent framework of Chiral Perturbation Theory (ChPT) [HI 
ini [ID]- It was noted in Ref. [S] that the radiative corrections to the ratio 1,2(1) /'^■n-^2{i) 
are uniquely determined to lowest non-trivial order in the chiral expansion and the final 
numerical results have been reported in Refs. [H] Il2]. So far, strong isospin breaking 
corrections oc (m„ — rrid) have been ignored in the extraction of |V^s/V^d|. In view of the 
increasing precision of lattice results on Pa'/Ptt, also these effects should be taken into 
account and it is the purpose of this letter to fill this gap. 

In Section [2] we perform our analysis of both electromagnetic and strong isospin viola- 
tion to rp£2(7) ill the unified framework of ChPT. In Section Isj we update the phenomeno- 
logical determination of |Kis/Kjd| and the global fit to l^irff^and \Vus\ [12] to reflect the 
new strong isospin breaking correction. We present our conclusions in Section |4j 

2 Isospin violation in tt ^ ^^^(t) and K ^^^(7) 

The fully photon-inclusive decay width 

Fp^^^^, =r(P^£z/,) + r(P^£z/,7) , P = 7r^,K^, (2.1) 
can be written in the form [H [9] 

rp,2(7) = ^^'4^''^^ Mp m\ [\ - ^Ew (1 + 5™) > K± = Vua , Vk^ = • 

(2.2) 
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Gp is the Fermi constant extracted from muon decay and Ftj± {Fx±) is the decay constant 
of the charged pion (kaon) in pure QCEQ Sew = 1-0232 is the universal electroweak short- 
distance enhancement factor [U El E] appearing in semileptonic decays. The long-distance 
electromagnetic corrections are described by |1] 

<^EM = ^ (Pime/Mp) + ^ In ^ - Cf + • • •) • (2-3) 

The one-loop function F{x) is related to the universal long-distance electromagnetic 
corrections (limit of point-like meson), its explicit form can be found in Ref. ^j. The di- 
mensionless constant C[ parametrizes electromagnetic effects due to hadronic structure. 
The dots refer to further structure dependent terms [9j arising at higher orders in the 
chiral expansion. 

The constants and Cf^ were calculated in chiral perturbation theory with virtual 
photons and leptons to order e^p^ |H]: 

= -E^M,) + - ^3 + 2 In ^ + lu ^ j , (2,4) 

C- = -£^(Ay + |(3 + 21ng+l„g) . (2.5) 

The parameter Z, appearing in the chiral Lagrangian of order e^p°, can be expressed 
in terms of the measured pion masses, the fine-structure constant and the pion decay 
constant: 

Z = ^ ^ ^ 0.8 . 2.6 

Both and depend on the same combination of electromagnetic low-energy cou- 
plings, 

E' = \+ 4vr^ i^-Kl + ^Kl + '^^K; + ^K^^ - \xi - AX^ + 4X3^ - Xf^') . (2.7) 
The Ki were defined in Ref. [15] and the Xi in Ref. |8j. Note that we have pulled out the 



short-distance enhancement factor in Eq. (2.2) keeping only the residu al lo ng-distance 



part Xq^^'^ jT6] in Eq. (2.7). The coupling constants entering in Eq. (2.7) have been 



estimated in Refs. [iTj and ^16j using large- Ai'c methods, giving 

= -2.4(5) , = -1.9(5) . (2.8) 

^In the case of the pion, the distinction between the charged and the neutral decay constant is a tiny 
effect of order (m^ — whereas Ff^±/Ff^o is of the order — r7i„ jl3j . For a discussion of some 

subtle points in the separation of QCD and electromagnetism, the reader is referred to Refs. [5], |14| . 
and [Tn|. Here we simply recall that the scale ambiguity inherent to the separation of QCD and QED 
effects in the meson decay constants in the chiral limit has been shown to be about 8 keV p!0^ , and thus 
negligible at the current level of precision. 
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The errors given here are based on naive power counting of unknown contributions arising 
at order e'^p^. 

For the ratio of decay widths we find 



r 



Mk± {1-mj/Ml^y 
\Vud\^ Fl± M^±{l-mi/Ml^f 



Jem J 



(2.9) 



where 
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Cf + CI 



F{m,/MK) - F{m,/M^) + 
-0.0069(17) . 



3-Z. M 



In 



K 



Ml 



(2.10) 



The important feature of the above result is that the electromagnetic low-energy coupling 
E'' cancels out in the ratio [8]. To lowest non-trivial order in the chiral expansion, the 
electromagnetic correction 5em is uniquely determined in terms of the fine-structure con- 
stant and masses, demonstrating the power of effective field theory methods. The quoted 
25% uncertainty is an estimate of corrections arising to higher order in the chiral expan- 



sion. Note also that in Eq. (2.9) the ratio Fj(±/F^± is the full QCD quantity including 



strong isospin breaking due to ruu 7^ m^. 

With experimental measurements of the 71^2(7) and K£2{-y) decay rates and the precise 
knowledge of the radiative corrections given in Eq. (2.10), Eq. (2.9) can be used to obtain 
the value of the ratio 



IK. I Fk± 



'Tk,^^^^M^± 



1/2 



mi/Ml 



IK 



ud 



1 - 5em/2) 



(2.11) 



It was suggested in Ref. [Ij to combine this result with lattice data on the ratio of the 
kaon and pion decay constant to determine |Ks/Kd|- Finally, using also experimental 
input for |Kd|, a value for |Ks| can be obtained in this way. 

As lattice calculations [IHl [IHl EQl EH |22l |23l El] are usually still performed in the 
limit of equal light quark masses (m„ = m^), it is convenient to rewrite Eq. (2.9) in the 
form 



<2(7) 



iKsP Fl Mk± (l-m2/M|±; 



r.,,(,, iKdP F^ M^^il-mj/MX 



1 + 5em + Si 



JSU{2) 



where F^^ and Fk denote the decay constants in the isospin limit and 
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5. 



SU{2) 



( Fk±F^ 
\F^±Fk 



- 1 



(2.12) 



(2.13) 



The isospin symmetry limit is defined by ruu = rrid, e = and an appropriate choice 
of the isospin symmetric meson masses M^, Mk- To the order we are working here, it 
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is sufficient to adopt the following convention: and Mk are related to the measured 
values of the masses by 



(2.14) 



as the leading order isospin violating contributions of strong and electromagnetic origin 
cancel in these expressions (see e.g. Ref |8]). For a fully consistent application of the 
isospin violating corrections calculated in this work, the lattice QCD results should be 
extrapolated to the isospin-limit "physical" meson masses and M^^. 

The chiral expansion of and up to the order p^, (m^ — mu)p^ is given by 



Fk± 



i^o 1 + 



2(47r)2Fo^ 
4 



8(47r)2Fo2 
8^3 £ 



L\{^,){Ml + 2Ml) + Ll{^^)Ml 

A/f 2 /\/f2 

2M,^ln^ + Miln^ 
LK/i)(M^2 + 2Mi) + L^5(/.)Mi 



^Ml In ^ + 6Mi In ^ + SM^^ In ^ 



2 A#2\ 



4(4^0W 



M. 



M^^ln-f - M^Mn-f - -(Mf - M^^) In 



M 



K 



(2.15) 



(2.16) 
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Fq is the pion decay constant in the limit of chiral SU(3) and M^^ and denote the 
isospin limits of the pion mass and the kaon mass, respectively, as discussed above. To 
lowest order, they can be expressed in terms of the quark masses by 



Ml = 2Bom , M| = Bo{m, + m) 



m 



(m„ + rrid) 



(2.17) 



where Bq is related to the quark condensate. To the same order, is related to Mk 
and M^ by the GMO relation 



M^ 



m 



-5o + - = -M^ - -Ml . 



(2.18) 



The parameter e of strong isospin breaking (which coincides with the lowest order 7r°-?7 
mixing angle) is given by 

' v^^'^-^^ (2.19) 



4 rris — 171 
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Taking e = in Eqs. (2.15) and (2.16) reproduces the results given in Ref. [TH] in the 
isospin hmit. The correction parameter (5su(2) can be read off as 



^SU{2) 



(2.20) 



Vs. 



2(4vr)2F2 



Ml 
^2 



Ml In — ^ 

/i2 



M 



V 



12- 



Expressing in terms of the isospin limit ratio Fk/F,^, (5su(2) can be recast into the 
following compact form: 



4 1 



Mi 



In^ 



(2.21) 



With the FLAG p6| averages of Nf = 2 + 1 lattice calculations, 



AR 



0.0116(13) , Fk/F^ = 1.193(6) 



(2.22) 
(2.23) 



we obtain 

^SU(2) = — 0.0043(5) (ll)higher orders , 

where we have estimated the uncertainty due to higher ordei[^ corrections in the chiral 
expansion to be at a level of 25%. We note that the strong isospin-breaking correction 
'5su(2) is of the same order of magnitude as the electromagnetic correction 6em in Eq. 
(2.10) and should not be neglected in the extraction of the ratio |Kis/Vud|. 



3 Update of and fit to |K 

We start from the basic formula 

1/2 1 o /nj2 



US 



\Vus\ 


Fk 


fVK,,,^,M^±\ 




F. 





mj/Ml^ (l - ^em/2 - 5su(2)/2) . (3.1) 



-1 



Putting together the results on isospin breaking obtained in the previous section and the 
experimental values for the leptonic widths of the pion '^■n^^^^-^ = 38.408(7) (/is) 
and the kaon Vk ^, = 51.25(16) (/is)~^ [11] we obtain: 



IK. 


Fk 


\Vud 


F. 



0.23922(25) x 



'p \ 1/2 



0.2763(5) . (3.2) 



An analysis of isospin violation beyond 0(p*) in ChPT [57] is beyond the scope of this work. 
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Figure 1: Graphical representation of current status of \Vud\, \Vus\ and corresponding 
CKM unitarity test. The horizontal band represents the constraint from from Ki^ decays, 
the thin vertical band the constraint from 0"^ — ?■ 0^ nuclear decays, the oblique band the 
constraint from K^2/'^^i2, and the ellipse is the la fit region. 



Finally, taking as reference value for Fk/F^^ the FLAG [2S] average of Nf = 2 + 1 



calculations (see Eq. (2.22)) we find for the ratio of the CKM elements: 

0.2316(12) . (3.3) 



\Vus\ 



K 



ud 



As discussed in Ref. [12], one can perform a fit to \Vud\ and \Vus\ using as input 
the values of \Vud\ = 0.97425(22) from superallowed nuclear f3 transitions [29], |Vus/Ktd| 



from Kg2/'n'e2 (Eq- (3.3)), and IV^^I from Kf^ decays. Measurements of Ki^ decays (BRs, 
form factors) [12], supplemented with theoretical calculations of electromagnetic and 
isospin corrections [30], EI] allow one to extract the product |Kts|/+(0) = 0.2163(5) [12] . 
Combining this with the Nj = 2 + 1 lattice result for the K ^ tt vector form factor 
/+(0) = 0.9599(34) (1^^) (14) [32] one obtains [26]: 

\Vus\iK,s) = 0.2255(5)exp(12)th ■ (3.4) 
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A fit to the above described input leads to 

iKdl = 0.97425(22) , |K.| = 0.2256(9) , (3.5) 

with x^/ndf = 0.012 and neghgible correlations between IV^^I and \Vus\- Fig. [l] provides 
a graphical representation of the various constraints in the \ Vus\ — I Kid I plane and the la 
fit region (in yellow). 

The fit values of \Vus\ and \ Vud\ (together with the neghgible contribution from \ Vub\ = 
0.00393(36) [llj) imply a stringent test of CKM unitarity: 

AcKM = iK^r + iK.r + iK^r - 1 = 0.0001(6) . (3.6) 

The remarkable agreement with the standard model prediction Ackm = allows one to 
set strong bounds on the effective scale of dimension-six operators that violate quark- 
lepton universality of charged current weak interactions [33j and thus parameterize new 
physics contributions to Ackm- The new physics effective scale is constrained to be A > 11 
TeV (90% C.L.), a lower bound to be compared to the effective scale of the standard model 
set by the Fermi constant, Asm = {2\/2Gf)^^^'^ — 174 GeV. This puts this charged- 
current quark lepton universality test at the same level as the precision electroweak tests 
from Z-pole measurements. 



4 Conclusions 

In this letter we have discussed in detail the effect of electromagnetic and strong isospin 
breaking relevant for the extraction of from the combination of the measured 

^K^2(j)/'^'^^2{-y) ratio of decay constants Fk/Ft^ calculated in lattice QCD in 

the isospin limit. We have performed our analysis in the unified framework of ChPT, 
providing analytic expressions and numerical estimates for the EM and strong isospin 



breaking corrections 6em and 5su(2) that appear in the ratio ^f^Tr 2( (^^^ -'^l- (2-12)) 



The main new result of our work is a compact analytic expression for the correction 



due to strong isospin breaking (see Eq. (2.21)). This effect, neglected in all previous 



extractions of is comparable in size to the electromagnetic correction: 

5em = -0.0069(17) , 5su(2) = -0.0043(12) . (4.1) 

Using the current experimental results for the decay widths and Fk/F^, = 1.193(6) [26] , 
we obtain |Vu<(/Vurf| = 0.2316(12), which is about 0.2% above previous determinations 
that ignored (5su(2)- 
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